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ABSTRACT 

The winds and radiation from massive stars clear out large cavities in the interstellar medium. These 
bubbles, as they have been called, impact their surrounding molecular clouds and may influence the 
formation of stars therein. Here we present JCMT observations of the J = 3 — 2 line of CO in 43 
bubbles identified with Spitzer Space Telescope observations. These spectroscopic data reveal the 
three-dimensional structure of the bubbles. In particular, we show that the cold gas lies in a ring, not 
a sphere, around the bubbles indicating that the parent molecular clouds are flattened with a typical 
thickness of a few parsecs. We also mapped 7 bubbles in the J = 4 — 3 line of HCO+ and find that 
the column densities inferred from the CO and HCO+ line intensities are below that necessary for 
"collect and collapse" models of induced star formation. We hypothesize that the fiattened molecular 
clouds are not greatly compressed by expanding shock fronts, which may hinder the formation of new 
stars. 

Subject headings: ISM: kinematics and dynamics — stars: formation — (ISM:) Hll regions 



1. INTRODUCTION 

Massive stars exhibit strong winds and radiation fields 
which clear out parsec-scale caviti es in the inter stellar 
medium (ISM) (Castor et al. 19751 : IWeaver et all [1977). 
At the boundaries of these cavities (or bubbles) lie the 
material displaced by the propagating wind or ioniza- 
tion front. Many have argued that such a dynamic pro- 
cess, which alters the physical environment of molecular 
clouds, might trigger the formation of a new generation 
of stars. Various mechanisms are reviewed bv lElmegreenl 
(11998,). 

Bubble surfaces are illuminated by the central stars 
and shine brightly in the infrared. A by-product of the 
Spitzer Galactic Legacy Infrared Mid-Plane S urvey Ex- 
traordinaire (glimpse; iBenjamin et al.l [20031) was the 
ident i ficatio n of nearly 600 bubbles (jChurchwell et al.l 
l2006l l2007l ). Because of Spitzer 's tremendous increase 
in sensitivity and resolution over previous infrared ob- 
servatories, the GLIMPSE bubble catalogs represent a 
dataset of interstellar bubbles unprecedented in size and 
detail. Simultaneously, sub-millimeter astronomy has 
enjoyed rapid technological improvements. The newly- 
commissioned HARP heterodyne receiver array on the 
James Clerk Maxwell Telescope (JCMT) provides the 
ability to efliciently map large portions of the sub- 
millimeter sky at high spatial and velocity resolution 
(|Smith et al.l I2OO8I ). This permits a followup study of 
the molecular gas in a sizable subset of the original bub- 
ble catalogs. 

In this paper we present HARP obser vations of the J = 
3 — 2 C O line toward 43 bubbles in the lChurchwell et al.l 
(|2006l ) catalog. The data reveal the three-dimensional 
structure of the bubbles and their interaction with the 
ambient molecular cloud. We also observed the J = 4 — 3 
HCO+ line toward 7 bubbles to determine the amount of 
dense gas in their surroundings. These data are analyzed 
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together with the glimpse catalogs, and with 20 cm ra- 
dio continuum images of the ionized gas, to quantify the 
degree of star formation in and around the bubbles. 

We find that the molecular gas tends to lie in rings 
and conclude that the bubbles quickly break out of their 
molecular surroundings. This, in turn, suggests that 
molecular clouds are oblate with one dimension signif- 
icantly thinner than the other two. We investigate the 
extent of cloud compression and star formation in these 
regions, and posit that triggered star formation may be 
hindered due to the central stars' rapid breakout from 
their molecular confines. 

In ^we describe our observations, and introduce the 
supplementary infrared and radio datasets we have uti- 
lized. We present and discuss our analysis of these data 
in ^ In 21 we interpret our findings in the context 
of current theories about molecular cloud structure and 
star formation. We conclude in ^ 

2. OBSERVATIONS 

HARP is a 4 X 4 heterodyne array on the JCM T de- 
signed for large scale mapping (iSmith et al.ll2008l ). Un- 
der typical atmospheric conditions on Mauna Kea, the 
HARP array can map 345 GHz line emission to a noise 
level of Trms = 0.3 K at a rate of 100 arcmin^ hour~^. 

We selected targets at Galactic longitude / > 10°, and 
with angular radii 2' < ^ < 10' from the Churchwell cat- 
alogs of 592 bubbles. These constraints are based on the 
practical considerations of being observable from Hawaii 
in the summer and being well resolved but not so large 
as to require a prohibitively long time to map. As the 
bubbles span a wide range in distance and size, there is 
little correlation between the angular and intrinsic size. 
Hence, our target list is a fairly uniform sampling of the 
original catalog. 

Our observations were carried out over ten nights in 
2008 July- August. Data were acquired in raster observ- 
ing mode, using a nearby reference position that was first 
verified to be free of emission. We fit and removed a 
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first order baseline to the spectral data away from the 
line, and then gridded and coadded the individual raster 
pointings. Finally, we binned the coadded spectral cubes 
to 0.2 km s~^, and smoothed the data with a Gaus- 
sian kernel. The final spatial resolution of each cube is 
IQ" . All data reduction was performed using the STAR- 
LINK ^ package. To convert the intensities to radiation 
temperatures, we adopted a main beam efficiency of 70%. 
In total, we observed 45 rings in the CO 3-2 transition, 
and 7 in IICO+ 4-3. The reduced data have typical rms 
values of Trms = 0.3K. 

To supplement our sub-mm data, we retrieved im- 
age and catalog data of these regi ons from the S pitzer 
GLIMPSE and VLA MAGPIS surveys (|Helfand et al.l [2006). 
The infrared data provide high resolution images of the 
bubbles, and help identify bubble features in the sub-mm 
data. Furthermore, we use the glimpse data to identify 
and study young stellar objects (YSOs) in these regions. 

The VLA MAGPIS survey is a multi-configuration 
VLA survey of cm -emission in the galactic plane 
(|Helfand et al.ll2006h . We use these data to identify Hll 
regions associated with our sample. The intensity of Hll 
regions provides an estimate of the ionizing photon fiux 
in the region, which in turn constrains the mass of the 
central star(s) driving bubble expansion. 

3. RESULTS 

3.1. Molecular rings 

CO is the most abundant molecule with a low-energy 
transition in the ISM. The J = 3 — 2 transition is an 
excellent tracer of the cool, 20-50K, moderately dense, 
^ 10^~^cm~^, gas around the bubbles. The emis- 
sion from this line traces the ambient molecular cloud 
around the bubble shells, and is particularly strong at 
their interface. HCO+ is moderately abundant but has 
strong emission in gas with densities > lO^cm"^. It 
is a useful signpost of star-forming and potentially pre- 
stellar gas. 

We detected CO emission toward 43 of the 45 bubbles 
that we observed and HCO+ in 6 of the 7 followup tar- 
gets. Images of the molecular emission in relation to the 
Spitzer 8 /im and, where available, 20 cm data are shown 
in Figure [TJ Figures for all 43 bubbles in our sample can 
be found in the electronic version of this article. 

There is generally a very good correspondence between 
the CO and Spitzer maps; because the interiors of bub- 
bles are well-evacuated, they are clearly defined in the 
sub-mm data cubes. There are occasional departures 
between the infrared and sub-mm morphologies. For ex- 
ample, filamentary structures can be seen in many of 
the sub-mm data, but do not emit prominently in the 
infrared. These structures extend both out of (e.g., bub- 
bles N37, N45, N90) and into (e.g., N29, N27, N44, N49) 
the bubbles themselves. 

The CO 3-2 line is readily excited by collisions in 
molecular clouds, and very easily becomes optically 
thick. Because of the CO's opacity, the jcmt data cubes 
show the surfaces (in position- velocity space) of molecu- 
lar clouds. The filamentary structures in these data illus- 
trate the inhomogeneity of molecular clouds on parsec- 
scales. These inhomogeneities are not as well probed at 

^ http:/ /starlink.jach. hawaii.edu/starlink 




Fig. 1. — Examples from our dataset. The left figures show the 
peak intensity of the CO and, when observed, HCO+ emission. 
The center figures show Spitzer 8/im emission in blue, and VLA 20 
cm emission in red. The rightmost figures show the first moment 
of the CO emission. The rings shown are N36 (top) and N45 (bot- 
tom). Equivalent figures for our entire sample can be found in the 
electronic version of this article. 

mid-infrared wavelengths. The emission at 8 /im is dom- 
inated by UV-exci ted PAH emission, and is expected to 
be op tically thin (|Churchwell et all l2QQ6l : I Wat son et al.l 
[2QQ8f) . The more uniform nature of bubble emission 
at these wavelengths most likely refiects the (approxi- 
mately) spherically symmetric UV radiation field pro- 
duced by the bubble-blowing stars. 

The most striking feature of the CO data is the rarity of 
emission towards the center of the bubbles. If bubbles are 
2-dimensional projections of spherical shells, we would 
expect to see the front and back faces of these shells at 
blueshifted and redshifted velocities. Of the 43 bubbles 
we have studied, we find no convincing examples of such 
a structure. 

Three different lines of evidence suggest that the 
molecular data are inconsistent with a spherical shell 
morphology. First is the contrast ratio between the bub- 
ble shells and bubble interiors. In the simplest case of 
optically thin emission, the integrated intensity scales 
linearly with the line-of-sight path length through the 
shell, as a function of impact parameter 6, 

' 2R[{1 + AR/R)^ - (h/Rf]^ -2R[l - {h/Rf]^ 

b<R 

^ " I 2R[{1 + AR/R)^ - (h/Rf] ^ R<b<R^AR 
,0 b>R^AR 

(1) 

where the bubble h as an inner radius R and thickness 
AR. Wat son et al.l (|2b08) have used this equation to 
describe the 8/im radial intensity profile for bubble N49. 

Compared to 8/im emission, CO provides a stronger 
test of this equation for two reasons. First, the added 
velocity information provides some ability to separate 
bubble emission from unassociated foreground and back- 
ground material. Second, the expected opacity of CO 
serves to decrease the edge-to-center contrast ratio, mak- 
ing the bubble fronts and backs easier to detect if present 
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Fig. 2. — Plot of the azimuthally averaged radial profile of four 
rings, compared to profiles for optically thin spherical shells. The 
data are inconsistent with these profiles, which have overly thin 
outer edges or low edge-to-center contrast ratios. This strongly 
disfavors the spherical shell model for bubble morphologies, and 
suggests instead that these bubbles lack fronts and backs. 



(the extreme profile is a completely opaque disk with no 
radial intensity variation). 

Figure [2] plots Equation [1] for three different values of 
/S.R/R. The azimuthally averaged radial profiles of four 
well defined bubbles are also plotted on this Figure. The 
Figure illustrates that the contrast ratio between the CO 
shells and interiors is problematically high. To reproduce 
such high contrasts, very thin shells must be invoked. 
However, the shell thicknesses are resolved in several of 
the larger objects in our sample, and rule out the very- 
thin shell hypothesis. Furthermore, as mentioned above, 
treating the CO emission as optically thin is an unre- 
alistic assumption, and its opacity makes the dearth of 
emission interior to bubble shells even harder to reconcile 
to a spherical shell morphology. 

The intensity of emission surrounding bubbles provides 
the second line of evidence against the spherical shell hy- 
pothesis. Most of the objects in our sample are found 
to lie interior to, or on the edge of, regions of extended 
molecular emission. If bubbles are embedded in molecu- 
lar material, then these clouds should contribute fore- 
ground and background emission inside bubble shells. 
However, the CO interiors of many bubbles are actually 
darker than their immediate exteriors (see, e.g., N21, 
N22, N29, N49, N74, N80). We note that even the 
8 jiia emission occasionally displays this signature (bub- 
bles from this study include N90, N94 and, to a lesser 
extent, N36 and N37). 

Finally, the expanding spherical shell hypothesis makes 
predictions about the velocity structure of bubbles which 
are not borne out by the CO data. An expanding spher- 
ical shell would manifest itself as a connected, coher- 
ent structure in a data cube. Specifically, as one steps 
through velocity space, a sphere would appear as a blue- 
shifted point (the bubble front) expanding into a ring 
(the midsection), and then contracting back to a red- 
shifted point (the back). A channel map for bubble 
N29 is displayed in Figure [3l The occasional filamen- 
tary structures seen interior to bubbles do not display 
the systematic velocity shifts as expected, but are in- 
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Fig. 3. — Example channel map for bubble N29. The radial 
velocity at each slice is listed in km s~^ in the upper right cor- 
ner. An ellipse is drawn in for reference, to mark the approximate 
boundary of the interior. Note that the bubble interior is largely 
devoid of emission. The few structures that do appear inside the 
shell (for example, at a radial velocity of 37.6 km s"-*^) are not 
red- or blueshifted, as would be expected if N29 was an expanding 
sphere. Note also that the extended emission to the north of the 
bubble does not extend into the interior of N29. This points to 
an absence of ambient molecular material in front of or behind the 
bubble. 



stead moving at roughly the same radial velocity as the 
edges. These filaments are either unaffiliated foreground 
or background structures (consistent with the observa- 
tion that they are absent from the Spitzer data), or they 
imply that the expansion of these regions has stalled. 

In light of this evidence, the most natural interpre- 
tation to draw from the observations is that the fronts 
and backs of bubbles are missing, and we are instead ob- 
serving only the bubble midsection. A corollary is that 
the molecular clouds in which bubbles are embedded are 
oblate with thicknesses of a few parsecs, comparable to 
typical ring diameters (Figure [4|). Expanding bubbles 
break out of the cloud along this axis and we see a ring of 
CO emission, approximately circular if the flattened axis 
is along our line-of-sight, otherwise elliptical. If viewed 
edge-on, the bubbles would not be identifiable as such ei- 
ther in our images or those from Spitzer, but might per- 
haps be classif ied as filanientary structures or Infrared 
Dark Clouds (Simon et al.ll2QQ6l : iJackson et al.l [2Q08). 

The relationship between the ring diameter and its par- 
ent cloud's thickness depends on the mechanism driv- 
ing the ring's expansion. If the stars powering the bub- 
ble produce stellar winds of sufficient strength, then the 
bubble expansion is driven by direct mome ntum transfer 
from the wind into the ambient material ([Castor et al.l 
I1975D . In this scenario, the ring's expansion would con- 
tinue after it breaches the flattened cloud, and the ob- 
served ring diameter places an upper limit on the thick- 
ness of the surrounding material. Alternatively, if the 
stars power an Hll region but do not have powerful 
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Fig. 5. — An illustration for how we defined each ring's bound- 
ary in the CO data cubes. The red and blue lines on the left 
trace, respectively, the intensity-thresholded and manually-defined 
boundaries for bubble N22; both boundaries are needed to eff'ec- 
tively define the inner and outer ring boundaries. The figure on 
the right is the intersection of these two regions, and defines which 
pixels we consider to be part of this bubble. 



Fig. 4. — A schematic of an interstellar ring, where a central 
group of O and B stars clear out a cavity in a fiattened molecular 
cloud. The top panel shows a face on view of an approximately- 
circular cavity of diameter a few parsecs. None, or very little, CO 
emission is seen towards the cavity, whereas the molecular cloud 
extends for tens of parsecs in the plane of the sky. We interpret our 
finding as indicating that clouds have line-of-sight depths of a few 
parsecs, and that bubbles rapidly break out of their host clouds 
(illustrated in the lower panel). In this scenario, the interaction 
region between the star and cloud is smaller than it would be than 
for a more homogeneous cloud geometry. 

winds, then a shock from the pressure-driven expansion 
of the Hll region s urrounds the stellar wind shock front 
(jFreyer et al ll2QQ3 ). Rings in this case trace the bound- 
ary of the Hll shock and, once the bubble breaches the 
cloud top and bottom, the pressure would drop and ex- 
pansion would stall. In this case, the ring diameter is 
roughly equal to the cloud thickness. As we discuss in 
Section 13. 3[ both classes of bubbles are present in this 
sample. 

The ring-like morphology of cold molecular gas en- 
circling spherically expanding structures appears to be 
commonplace: it is present not only in the data pre- 
sented here, but also in many other well studied examples 
of large Hll regions ( Willi ams et a l. 1995; Koenig et al. 
l2QQ8l : IZavagno et al.ll2QQ6l : iDeharveng et al.ll2QQ 9). How- 
ever, the implications of this observational signature have 
not generally been commented upon. Furthermore, these 
bubbles are smaller than typical optically visible Hll re- 
gions, and thus place stronger constraints on their parent 
clouds' thicknesses. 

3.2. Ring properties 

While the molecular rings can often be easily identified 
by eye in a CO data cube, it is harder to rigorously define 
where the bubble interaction zone ends and the ambient 
molecular cloud begins. Here we describe our process of 
isolating bubble features in the data. 

We first cropped the data in velocity to include only 
those channels with matching features in the Spitzer im- 
ages. In most cases, the inner ring edge stands out in the 
integrated maps. Thus, the inner boundary is defined by 
an integrated intensity threshold, chosen manually for 
each ring. While some outer ring boundaries stand out 
against dark backgrounds in the integrated maps, rings 
on the edges of extended cloud structures do not. As 



a compromise, we define each ring's outer boundary as 
the intersection of the intensity threshold and a manual 
boundary. An example of this process is illustrated in 
Figure O the red line on the left traces the boundary 
defined by hand, while the blue line traces the intensity 
threshold; the figure on the right shows the intersection 
of these two regions, and serves to define which pixels 
are associated with this ring. 

The ring properties for our 43 CO detections are listed 
in Table [H Near kinematic dista nces are calculated us - 
ing the galactic rotation model of lBrand fc Blitzl (1993 ), 
for Rq = 8.8 kpc and Vq = 275 km/s (jReid et al.ll2QQ9[ ). 
Kinematic distance uncertainties are dominated b y pecu - 
liar motions. Based on the data from Reid et al. (l2009l), 
we estimate this uncertainty to be about 30%, and use 
t his value for the d i stance errors reported in Table [H 

IChurchwell et al.l (|2QQ6l ) make the argument that bub- 
bles are more likely located at their near kinematic dis- 
tances, since interstellar extinction and confusion with 
other diffuse emission structures would tend to obscure 
objects on the far side of the Galactic disk. Following 
this assumption, we convert the Spitzer angular shell 
radii and thicknesses to linear measurements, and include 
them in Table [H 

Seventeen of the bubbles i n our CO survey are associ- 
ated with the Hll regions in iPaladini et all (|2QQ3[ ). Our 
independent velocity measurements agree with the re- 
combination line velocities in that data set, strength- 
ening the argument - made previously on statistical 
grounds - that features seen in the infrared and sub- 
millimeter are physically associated with these ionized 
regions. Furthermore, we provide 26 new kinematic dis- 
tance measurements. We note that submillimeter data 
have the advantage of providing distance estimates to 
rings without detected Hll regions. 



3.3. Constraints on the central sources 

The stars that create bubbles also ionize the material 
in their immediate vicinity, which glows via recombina- 
tion and free-free processes. The VLA MAGPIS survey is 
sensitive to this emission, and can constrain the types of 
stars responsible for the ionization. The MAGPIS survey 
overlaps 40 bubbles in our sample. 
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TABLE 1 
Bubble properties 



Bubble^ 


I 


b 






C^near 










(°) 


(°) 


/I 1 \ 
(kms ^) 


/I 1 \ 

(kms ^) 


(kpc) 


(kpc) 


(pc) 


(pc) 


N5 


12.46 


-1.12 


39.5 


2.4 


3.7 


1.1 


3.49 


0.71 


NIO 


13.20 


0.06 


50.2 


4.1 


4.1 


1.2 


1.57 


0.42 


N14 


13.99 


-0.13 


40.3 


2.8 


3.5 


1.1 


1.24 


0.39 


N15 


15.01 


-0.61 


19.7 


5.2 


1.9 


0.6 


0.83 


0.25 


N16 


14.98 


0.06 


25.2 


3.0 


2.4 


0.7 


1.73 


0.44 


N20 


17.92 


-0.67 


42.1 


3.3 


3.2 


0.9 


1.04 


0.25 


N21 


18.19 


-0.39 


48.1 


3.1 


3.4 


1.0 


2.14 


0.49 


N22 


18.26 


-0.30 


51.3 


2.1 


3.6 


1.1 


1.77 


0.45 


N27 


19.81 


0.02 


57.6 


2.6 


3.8 


1.1 


1.19 


0.24 


N29 


23.06 


0.56 


37.9 


2.1 


2.5 


0.8 


1.98 


0.59 


N30 


23.11 


0.58 


38.7 


2.6 


2.6 


0.8 


0.81 


0.26 


N34 


24.30 


-0.17 


101.0 


1.8 


5.2 


1.5 


1.65 


0.41 


N35 


24.50 


0.24 


115.0 


6.1 


5.7 


1.7 


5.41 


1.54 


N36 


24.84 


0.11 


109.0 


5.2 


5.5 


1.6 


4.51 


1.39 


N37 


25.30 


0.30 


40.8 


1.6 


2.6 


0.8 


1.34 


0.37 


N39 


25.36 


-0.15 


63.8 


3.8 


3.7 


1.1 


2.14 


0.70 


N40 


25.37 


-0.36 


54.8 


3.1 


3.3 


1.0 


1.18 


0.36 


N44 


26.83 


u.oo 


ox . X 


1.8 


4.3 


1.3 


1.41 


0.30 


N45 


26.99 


-0.05 


83.2 


3.4 


4.4 


1.3 


1.88 


0.68 


N46 


27.31 


-0.12 


92.6 


2.8 


4.8 


1.4 


1.93 


0.67 


N47 


28.03 


-0.03 


99.2 


4.4 


5.1 


1.5 


3.32 


0.76 


N49 


28.83 


-0.23 


87.5 


3.1 


4.6 


1.4 


1.77 


0.43 


N50 


28.99 


0.10 


70.2 


2.4 


3.8 


1.1 


1.88 


0.45 


N51 


29.15 


-0.26 


94.7 


3.0 


4.9 


1.5 


2.79 


0.60 


N52 


30.74 


-0.02 


94.1 


9.6 


4.8 


1.5 


2.85 


0.87 


N53 


31.16 


-0.14 


41.6 


2.6 


2.4 


0.7 


0.59 


0.14 


N54 


31.16 


0.30 


102.0 


4.8 


5.2 


1.6 


2.69 


0.57 


N56 


32.58 


0.00 


100.0 


2.8 


5.2 


1.6 


1.69 


0.42 


N61 


34.15 


0.14 


57.3 


3.5 


3.1 


0.9 


3.15 


0.51 


N62 


34.33 


0.21 


56.2 


4.2 


3.1 


0.9 


1.37 


0.30 


N65 


35.02 


0.33 


52.5 


2.6 


2.9 


0.9 


1.81 


0.46 


N74 


38.93 


-0.40 


40.4 


2.6 


2.3 


0.7 


0.92 


0.31 


N77 


40.42 


-0.05 


68.8 


2.0 


3.9 


1.2 


1.46 


0.36 


N79 


41.52 


0.03 


14.1 


3.1 


0.7 


0.2 


0.26 


0.06 


N80 


41.93 


0.03 


17.7 


2.8 


0.9 


0.3 


0.48 


0.11 


N82 


42.11 


-0.62 


66.5 


1.8 


3.8 


1.1 


1.86 


0.49 


N84 


42.83 


-0.16 


61.2 


1.2 


3.5 


1.1 


1.18 


0.31 


N90 


43.76 


0.09 


69.4 


2.3 


4.1 


1.2 


2.08 


0.50 


N92 


44.33 


-0.83 


62.0 


2.6 


3.7 


1.1 


2.09 


0.57 


N94 


44.81 


-0.48 


47.0 


2.1 


2.7 


0.8 


2.95 


0.91 


N120 


55.27 


0.25 


29.3 


2.6 


1.9 


0.6 


0.76 


0.17 


N130 


62.37 


-0.54 


0.4 


0.3 


8.4 


2.5 


2.57 


0.64 


N133 


63.15 


0.44 


21.1 


3.1 


1.6 


0.5 


0.76 


0.18 



^from Churchwell et al. (2006) 
^assuming near kinematic distance 



While free- free emission is clearly evident in the major- 
ity of our sample (see Figure [1]), there are several bub- 
bles for which radio emission is faint or hidden in the 
background. For both bright and faint Hll regions, we 
derive a flux estimate in the following way. We use aper- 
ture photometry to compute the integrated signal within 
each bubble, corrected for any nonzero background level. 
We refer to this signal measurement as because it is 
derived from the data. Because we find that the back- 
ground noise in the MAGPIS maps is well described by 
a Gaussian distribution, the error on a spatially inte- 
grated flux measurement is also distributed as a Gaus- 
sian, with standard deviation a = aoy^Tt/TtB' Here, ctq 
is the standard deviation of the background noise, ft is 
the solid angle subtended by the object of interest, and 
is the solid angle subtended by the antenna beam. 
Using Bayes' Theorem, the posterior distribution for a 



source's true integrated flux F is then given by 

P{D\F)P{F) 



P{F\D) = 



C^P{D\F)P{F)dF 



(2) 



Where P{D\F) is the likelihood measuring D when the 
true flux is F; it is simply a Gaussian centered on F, 
with a width of a, evaluated at D. The function P{F) is 
the prior probability distribution for F before any data 
were taken; we deflne this function to be zero for F < 
(since no real source can produce a negative flux), and 
constant otherwise. The actual value of the constant is 
irrelevant, since it cancels with the denominator. 

The beneflt of deriving flux estimates from Equation [2] 
is twofold. First, the prior P{F) naturally incorporates 
the the physical restriction that the true flux be non- 
negative. Second, it makes no distinction between de- 
tections and non-detections, and is equally valid in both 
situations. For each source in the MAGPIS data, we re- 
port the median of P{F\D) as the 6 cm flux estimate, 
since the true flux lies above or below this value with 
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equal probability. For strong detections, this value con- 
verges to D. For non-detections, the data provide upper 
limit constraints more effectively than lower limits, and 
the median of P{F\D) is comparable to a. 

These fluxes are given in Table [2] for the 4 bubb les cov- 
ered by the MAGPIS survey. Condon fc RansomI ([2006) 
derive a relationship between this flux and the ionizing 
photon rate: 

( - V V ^ (3) 

VGHz; VlO^o WHz-V 
this rate is provided in Tabled along with the number of 
09.5 stars needed to produce such a flux (Ma rtins et al.l 
[2005). These flux estimates are lower limits, since any 
ionizing photons that are absorbed by dust or escape 
out of the bu bble fronts or backs a re unaccounted for in 
Equation [3l IWatson et all (|2008[ ) identified the likely 
stars powering three bubbles, and estimated that the 
MAGPIS estimate of is about a factor of two lower 
than the expected ionizing photon flux from these stars. 
Thus, the values listed in Table [2] should be considered 
factor of ~ 2 underestimates of the true ionizing lumi- 
nosity. 

The Hll regions in this sample are substantially weaker 
than classical Hll regions. For comparison, we include 
data on the Orion Nebula (considered faint by classi- 
cal Hll st andards) and NGC 3603 (a giant Hll region, 
iKennicut t 1984). With a few exceptions, most of the 
free-free emission in this sample is fainter then even the 
Orion Hll region. Most of the bubbles are powered by 
early B stars, and can be viewed as lower-mass equiva- 
lents to classical Hll regions. This confirms Churchwell et 
al's (2006) original interpretation that bubbles surround 
late O and early B stars. However, we note that many 
of these objects do have some detectable level of free-free 
emission, and as such may be driven by the expansion of 
the Hll region as opposed to stellar winds (|Frever et al.l 
te.). 

3.4. Properties of the dense gas 

We mapped the HCO+ 4-3 line toward 7 of the 45 
bubbles to constrain the density of the gas in the bubbles. 
As the frequency of this line, 357 GHz, is very similar to 
that of CO 3-2, maps of each line have similar resolutions 
(16^') and can be directly compared. 

HCO+ was detected (Tpeak > 1 K) in 6 of the 7 bub- 
bles. Its detection indicates high density gas nu^ > 
5 X 10^ cm~^ (lEvans 1999). However, it is only seen 
in a single, generally small, part of the molecular ring 
where the CO emission was bright, Tpeak ^ 20 K. In sev- 
eral cases these dense regions coincide with CO outflows 
and, in one case (N74), a YSO overdensity, indicative of 
recent star formation (see also §3.5.ip . Table [3] lists the 
peak temperatures of the CO 3-2 and HCO+ 4-3 tran- 
sitions. Except for bubble 45, which was not detected, 
there are two entries for each bubble showing the values 
where the integrated intensity of HCO+ is greater than 
IKkms"^ and averaged over the entire ring as defined 
in go 

We compared these observations with radiative trans- 
fer models via the Large Velocity Gradient formalism 



TABLE 2 

H II REGION PROPERTIES 



Bubble 


20 cm Flux 

(Jy) 


log ^Ly 
(log s-i) 


JN (09. 5j 


N5 




47.13 


4 


NIO 




48.73 


15 


N14 


9 41 p_i_nn 


48.36 


6.4 


N15 




45.75 


.02 


N16 


1 7'^p_n'^ 

-L . ( Oc UO 


44.89 


.002 


N20 




46.87 


.20 


N21 




48.43 


7.3 


N22 




48.55 


9.8 


N27 


c; n4p— 09 


46.76 


.16 


N29 


4 Q4p_Lnn 


48.33 


5.8 


N30 




48.27 


5.1 


N34 




47.70 


1.4 


N35 




49.28 


53 


N36 




49.34 


60 


N37 


7 04e— 01 


47.57 


1.0 


N39 


1 fiOp-k01 


49.24 


47 


N40 




47.65 


1.2 


N44 




47.10 


4 


N45 


1 49p_i_nn 


48.33 


5.9 


N46 


1 OQp-kOO 


48.29 


5.4 


N47 




48.56 


10 


N49 




48.21 


4.5 


N50 


6 90p— 01 


47.85 


1.9 


N51 


9 qsp_n9 


46.75 


.15 


N52 


7 X4g_^01 


50.11 


350 


N53 


4.71e— 01 


47.33 


.58 


N54 


q 9qp_ni 


48.29 


5.4 


N56 


X.46e— 01 


47.49 


.85 


N61 


3.62e-01 


47.43 


.75 


N62 


5.87e-02 


46.64 


.12 


N65 


5,76e-02 


46.58 


.10 


N74 


2.56e-02 


46.03 


.03 


N77 


1.68e-02 


46.30 


.06 


N79 


3.94e-01 


46.18 


.04 


N80 


2.54e-01 


46.21 


.04 


N82 


1.39e+00 


48.20 


4.3 


N84 


5.20e-03 


45.70 


.01 


N88^ 


1.26e+01 






N90 


1.61e-01 


47.33 


.58 


N94 


1.19e-02 


45.83 


.02 


Orion^ 




48.85 


20 


NGC 3603^ 




51.05 


3100 



^No distance information (not detected in CO data) 
^Data from lKennicuttI (119841) 



using the program LVG in the MIRIAD^ package. For 
a given gas kinetic temperature, this program calculates 
the line radiation temperature of a specified transition as 
a function of the molecular column density (normalized 
by the linewidth) and H2 volume density. We created 
model grids of CO 3-2 and HCO+ 4-3 over the density 
range n(H2) = 10^ — 10^"^ cm~^ and normalized column 
density A/'(CO)/Av = 10^^ - 10^^cm-2(km s-^)-^ and 
7V(HC0+)/Av = 10^^ - 10i^cm-2(km s-^)-\ At these 
high densities, the 3-2 level of CO is in collisional equi- 
librium and the line strength is essentially independent 
of n(H2). It is mainly dependent on the normalized col- 
umn density, Nqo/Av. As the 4-3 level of HCO+ has a 
higher critical density, however, its line strength can be 
used as a diagnostic of n(H2) in this range. 

The models require kinetic temperatures Tkin > 40K 
to match the brighter bubbles (especially N29 and N82). 
These high values are not surprising given that the gas 

^ http:/ /carma. astro. umd.edu/miriad/ 
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lies on the edge of a wind-swept ionized region around 
massive stars. We further find that the inferred values of 
column and volume density are relatively insensitive to 
Tkin in the range 40 — 60 K. 

We plot the observed CO and HCO+ line temperatures 
on the LVG model grid for Tkin = 50K in Figured The 
filled circles show the regions of strong HCO+ emission 
and the open circles show the averages over the rest of the 
bubble. The two are clearly offset indicating, not surpris- 
ingly, both lower column and volume densities away from 
the HCO+ emission. We multiply the inferred values of 
N{CO)/Av by the observed CO linewidth and divide by 
an abundance [C0]/[H2] = 10~^ to derive the column 
density of the gas in the bubbles. The inferred volume 
density depends on the relative abundance of HCO^ to 
CO. Figure E] plots the models for [HCO+]/[CO] = lO""^ 
but a lower abundance would shift the HCO+ (dashed) 
lines up relative to the CO (solid) lines. In this case, 
the data points would move to the right implying higher 
H2 volume densities. The CO abundance is fairly well 
determined since it self-shields against dissociating radi- 
ation and is not depleted in the warm gas that we are 
observing. The HCO+ abundance is less certain, how- 
ever, especially as the bubbles are effectively large photo- 
dissociation regions. Nevertheless, the relative column 
densities of HCO+ (e.g., comparing two different bub- 
bles, or comparing a region of bright HCO^ emission to 
the average over the shell) are still well characterized. 
The values inferred from Figure [6] are listed in Table [3] 
and plotted in Figure [71 

We find column densities N(R2) = 5 — 10 x 10^^ cm~^ 
averaged over the regions with HCO+ emission. These 
are higher than the average over the rest of the bubble 
by a factor ~ 3 — 7. The volume densities, n(H2) range 
from 5 X 10^ to 5 x 10^ cm~^ and are enhanced by about 
an order of magnitude relative to the average in the rest 
of the bubble. 

In exp loring the conditions ne cessary for gravitational 
collapse. I Whitworth et al] (|l994l ) show that, in a wide va- 
riety of scenarios (cloud-cloud collisions, supernova rem- 
nant expansion, wind-blown bubbles, and Hll regions), 
hydrogen column densities of 4 — 8 x 10^^ cm~^ are re- 
quired for layers of molecular gas to collapse. While some 
of the dense knots of bright HCO+ in our sample may be 
unstable to gravitational collapse, the shells on average 
are an order of magnitude too tenuous to collapse due to 
gravitational instabilities. 

3.5. Star Formation 

As they sweep up and potentially compress ambient 
material, the expanding shock fronts from interstellar 
rings may drive new epochs of star formation. Here we 
examine the evidence for triggered star formation pro- 
vided by outfiows and overdensities of Young Stellar Ob- 
jects (YSOs). 

3.5.1. Young Stellar Objects 

With photometry alone, it is difficult to determine 
whether a given object is embedded in a bubble shell 
or merely aligned along the line of sight. Because of this, 
associated YSOs must be identified in a statistical sense 
by looking for YSO overdensities due to regions of en- 
hanced star formation. 



16.5 r 



12.5 



5.0 



5.5 



6.0 



6.5 

n(H2) 



7.0 



7.5 



Fig. 6. — LVG model plots of CO and HCO+ peak temperature 
as a function of column and number density. A kinetic temperature 
of Tk = 50K is assumed, along with abundances [CO]/[H2] = 10~^ 
and [HCO+]/[H2] = 10~^. The solid lines trace contours of peak 
CO line intensity, while the dashed lines trace contours of peak 
HCO"^ intensity. The filled circles correspond to regions of bright 
HCO"^ emission, whereas the open circles correspond to the average 
over the entire bubble. The open circles are upper limits to the 
HCO+ peak temperature. Figure [Jl gives which colors correspond 
to which bubbles. 




16 29 36 45 65 74 82 
Bubble 

Fig. 7. — The number and column densities for H2, inferred from 
LVG modeling. As in Figure (6] the filled circles are averages over 
regions of bright HCO+ emission, whereas open circles are averages 
over the entire shell. The open circles in the lower plot are upper 
limits, since HCO"^ is not detected in the average shell spectrum, 
and CO emission does not constrain the density. 



The GLIMPSE Highly Rehable point source catalogs 
(GLMIC) provide IRAC colors of point sources toward 
the bubbles. We have used this catalog to identify YSOs 
towards the bubbles in our sample. To identify likely 
Y SOs, we have use d the series of color cuts described 
in lGutermuth et"al] (|2QQ8l ). These cuts attempt to dis- 
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TABLE 3 
Properties of the shell emission. 





-rCGglOn 


rppeak 

(K) 


lAVcO 

(km s ^) 


rppeak 
HCO + 

(K) 


(km s -^j 


log nH2 
(log cm "^j 


log IN ynz J 
(log cm 


N16 


c 


5.98 


6.04 






5.81 


20.28 


N16 


H 


12.7 


6.68 


1.20 


2.96 


6.52 


20.69 


N29 


C 


14.0 


4.23 






5.45 


20.53 


N29 


H 


27.7 


3.22 


1.11 


2.51 


5.86 


20.83 


N36 


C 


7.98 


10.0 






6.11 


20.60 


N36 


H 


16.8 


10.0 


1.59 


8.57 


6.92 


21.06 


N45 


C 


4.84 


7.08 






5.66 


20.25 


N65 


C 


10.7 


6.13 






5.96 


20.51 


N65 


H 


16.5 


8.63 


3.24 


6.98 


7.48 


21.05 


N74 


C 


15.3 


5.00 






6.11 


20.66 


N74 


H 


20.9 


5.92 


2.39 


4.32 


6.67 


20.94 


N82 


C 


11.2 


3.77 






5.56 


20.34 


N82 


H 


24.4 


2.33 


0.65 


2.75 


6.01 


20.64 



indicates that the flux was averaged over the region of bright HCO"^ emission. C indicates 
that the flux was averaged over the entire shell, deflned by bright CO emission as outlined in 



go 





TABLE 4 

CO OUTFLOWS 



Outflow Bubble 







n 


n 


1 


N5 


12.405 


-1.132 


2 


N5 


12.427 


-1.112 


3 


NIO 


13.177 


0.061 


4 


N20 


17.932 


-0.664 



(k 



Fig. 8. — The two bubbles in our sample with YSO overdensities 
on the shell periphery. The positions of sources flagged as likely 
YSOs are shown as points. N74 has an extended cluster of YSOs at 
its southeast boundary, while N133 has a tight grouping of YSOs 
on its southern shell . Th e contours trace the bubble boundaries as 
discussed in Section [22] 

cern between bona- fide YSOs, star forming galaxies and 
AGN, and non-stellar PAH emission. We apply the cuts 
on de-reddened magnitudes; reddening is estimated us- 
ing the 2MASS-der i ved a ll sky extinction maps from 
iRowles fc FroebrichI (|2QQ9l ). We then examined each im- 
age for clusters of likely YSOs which coincide with bubble 
features. 

We find two clear examples of YSO overdensities, 
which we present in Figure [H The majority of bubbles 
have a few likely YSOs on their shells, but at a level com- 
parable to the field YSO surface density. Six bubbles 
shells (N15, N21, N27, N34, N39, N44) do not overlap 
any likely YSOs. 

In total, this method yields 109 YSOs which overlap a 
bubble shell as defined in §3.21 This represents a moder- 
ate enhancement over that expected from the mean YSO 
surface density towards these regions (which predicts 55 
chance overlaps). 

3.5.2. Outflows 

Young protostars have strong outflows that are often 
apparent as broad CO line wings. We searched the mo- 
ment maps for such features in the molecular rings and 
found 12 candidate objects. Their position and velocity 



'^low '^high 

ms"-*-) (kms"-*-) 



25 47 

35 45 

40 50 

34 46 



5 


N34 


24.295 


-0.153 


85 


115 


6 


N36 


24.850 


0.084 


103 


116 


7 


N36 


24.788 


0.082 


90 


135 


8 


N37 


25.285 


0.266 


30 


40 


9 


N40 


25.389 


-0.368 


30 


54 


10 


N49 


28.830 


-0.254 


80 


100 


11 


N65 


35.020 


0.348 


45 


70 


12 


N74 


38.934 


-0.362 


18 


56 



range are listed in Table [H Since this is an eye-based 
identification, the list is biased towards outflows that are 
more easily seen on top of the extended shell emission 
(e.g., strong outflows oriented towards our line of sight). 
At the ~ 0.25 pc resolution of our data, we did not resolve 
any clear bipolar structures. The outflow associated with 
N74 also coincides with the YSO overdensity there. 

3.5.3. Extended Green Objects 

'Cvgano wski et aP (|2008l ) recently identifled several 
hundred Extended Green Objects (or EGOs) in the 
GLIMPSE images. These objects are resolved structures 
which emit brightly at 4.5/im, and are believed to trace 
shocked molecular gas due to outflows from massive 
young stars. Three of these objects overlap the molec- 
ular shells of bubbles N39, N49, and N65, and two of 
these three (those associated with N49 and N65) are also 
detected as outflows in CO. 

4. DISCUSSION AND COMPARISON WITH 
PREVIOUS WORK 

4.1. Rings and Molecular Clouds 

The ring morphology of the CO emission implies that 
the molecular clouds in which bubbles form are flattened, 
with thicknesses not greater than the bubble sizes of 
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a few pc. Inferring the three dimensional structure of 
molecular clouds is a difficult end eavor, and previou s ef- 
forts have yielded mixed results. iKim et alT (|2QQ8l ) use 
supernova light echoes to probe the structure of the ISM 
behind Cas A. They observe numerous filamentary and 
sheetlike structures, with a low ffiling factor of ~0.4%. 
However, their technique filters out spatial structures 
larger than ~2 pc, and their data probe only the sub- 
structures in a potentially more extended region. 

Statistically analyzing a set of 10^ clouds. iKerton et al.l 
(2003) compare the observed distribution of molecular 
cloud axis ratios to those expected from triaxial ellip- 
soids, and conclude - in conflict with our result - that 
molecular clouds are prolate. The results of this tech- 
nique, however, are sensitive to the choice of model fig- 
ures (triaxial ellipsoids in Kerton et al.'s analysis). For 
example, the same technique has been applied to studies 
of molecular cores. Using two differe nt classe s of shapes 
(axisymmetric vs. tr iaxial ellipsoids). iRvdenI (pTOQG) and 
Uones & Basu (2002) come to opposite conclusions about 
whether molecular cores are prolate or oblate. The flat- 
tened nature of molecular clouds implied by our data is 
independent of these kinds of model assumptions, and is 
a direct inference from the observed lack of bubble fronts 
and backs. 

Our CO maps were limited to the bubbles and their 
environs. However, in many cases, we see extended emis- 
sion to the edges of our maps. In general, the bubbles 
lie in moderate-sized molecular clouds with dimensions 
greater than 10 pc in the plane of the ring. That is, the 
clouds have aspect ratios of at least a few, and perhaps 
as much as 10. 

Sheet-like ISM morphologies may naturally arise 
in molecular clouds created by the turbulent colli- 
sion and compression of warm neutral flows in the 
Galaxy (see, e.g., numerical work by Heitsch et al. 2005, 
IVazquez-Semadeni et al.l I2006L and references therein). 
The pervasive oblateness of the bubble-hosting clouds 
may be relevant to studies of molecular cloud creation 
and, by extension, star form ation (see recent reviews by 
[McK ee & Ostriker 20(51 and ^nnecker & Yorke 200Z). 

4.2. Rings and Triggered Star Formation 

While studies of individual "triggered" or "induced" 
star formation regions abound in the literature (recent 
papers on triggered st ar formation around bubbles and 
Hll regions in clude D eharveng et al.l 12 009. Snider et al. 

l2009LlKoeni7et al 20081 . and lPomarSs et al. 2009 ). there 

is less consensus on the extent to which the varied trigger- 
ing mechanisms contribute to the net stellar production 
in our Galaxy. 

The observation that molecular gas lies in a thin ring 
around bubble-blowing stars may well impact the abil- 
ity for these bubbles to induce new generations of star 
formation. In such a configuration, Hll regions likely 
lose pressure - and their ability to compress molecular 
gas - upon breaking out of the cloud. Similarly, momen- 
tum driven winds encounter a much smaller solid angle of 
molecular gas to sweep up. We have cataloged four trac- 
ers of future or ongoing star formation activity towards 
these regions: dense gas revealed by HCO+ emission, CO 
outflows, over densities of Young Stellar Objects, and co- 
incident Extended Green Objects. Here we examine to 
what extent these data constrain estimates of triggered 



star formation activity. 

The low detected fraction of dense HCO+ gas suggests 
that bubble shells are unstable to gravitational collapse 
only in small regions; the shells as a whole do not seem 
to be undergoing gravitational collapse. 

Because bubbles are powered by early B and late O 
stars (M > lOM©), any self-sustaining triggered star 
formation in these regions must produce stars massive 
enough to blo w their own b ubbles. Assuming a typi- 
cal IMF (e.g. Muench et al. 2002), an equivalent state- 
ment is that bubble-induced star formation must produce 
> 250 new stars in order to form a star massive enough 
drive a secondary bubble. 

The completeness of the YSO sample extracted from 
the GLIMPSE data is lim ited by t he bright nebular back- 
ground. Studying M17, iPovich et al. (2009) found that 
the GLIMPSE Point Source Archive recovers essentially 
all YSOs above ~ 3Mq towards these regions, while the 
completeness falls to - 10% at M = 2Mq, and - 1% 
at M = IMq. Our sample is less complete, as we have 
restricted our focus to the smaller Highly Reliable Cata- 
log, and most of our bubbles are 1.5-3 times more distant 
than M17. 

Two bubbles show suggestive YSO over densities, con- 
sisting of roughly 10 and 20 associated objects. If we 
assume, for the sake of simplicity, that no objects less 
massive than 3 M© are detected, then each flagged 
Y SO represents about 1 in 50 total stars formed via 
a iMuench et aD (|2002r ) IMF. Hence, it is possible that 
these overdensities trace clusters of > 500 — 1000 mem- 
bers, which would be sufficiently large to create new gen- 
erations of stellar bubbles. 

The handful of EGOs and CO outflows found towards 
these bubbles further suggests some degree of ongoing 
star formation in these regions. Unfortunately, it is dif- 
ficult to draw conclusions about the extent of star for- 
mation suggested by these objects. Molecular outflows 
are a feature of both high and low mass star formation, 
and do not constrain the number or mass of stars being 
formed. While EGOs are believed to trace specifically 
massive star formation, there are too few detections to 
draw any statistically meaningful conclusion; the overlap 
of EGOs and bubble shells is higher than that expected 
from chance alignments (3 alignments compared to .35 
expected), but the Poisson uncertainties on this result 
are very large. 

In all, up to 12 of the 43 bubbles in our sample may 
be undergoing c urrent star formati on as suggested by the 
above evidence. IChurchwell et al] ((2007) also looked for 
triggered star formation around bubbles (by looking for 
secondary bubbles, and YSO coincidences), and found 
evidence of such around 12% of the objects. 

Identifying triggered YSOs inside bubble shells is com- 
plicated by the fact that these objects are smaller 
and more distant than more h eavily studied regions 
(|Koenig et al.ll2008l : [Povich et al.ll2009f ). Furthermore, it 
is difficult to conffim that a given YSO was formed as the 
result of a triggering event as opposed to forming inde- 
pendently nearby. The soon-to-be-operational SCUBA2 
instrument on the JCMT will circumvent the limitations 
inherent to this method by providing sensitive measure- 
ments of the cold gas and dust column density towards 
these bubbles. These data will better determine whether 
or not these rings are gravitationally unstable, and hence 



10 



Beaumont & Williams 



whether the physical conditions of interstellar bubbles 
are conducive to star formation. 

5. CONCLUSIONS 

We have presented new CO 3-2 maps of 43 Spitzer 
identified bubbles in the Galactic plane. These maps 
point towards a new interpretation for bubble morphol- 
ogy - namely, that the shock fronts driven by massive 
stars tend to create rings instead of spherical shells. We 
suggest that this morphology naturally arises because the 



host molecular clouds are oblate, even sheet-like, with 
thicknesses of a few pc, and widths > 10 pc. Such a mor- 
phology implies that expanding shock fronts are poorly 
bound by molecular gas, and that cloud compression by 
these shocks may be limited. This idea is borne out by 
observations in HCO+, which show dense gas to be con- 
fined to small regions of bubble shells. 

We thank Ed Churchwell for stimulating discussions 
and the NSF for support through grant AST-0808144. 
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